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ABSTRACT 

This is the first in a series of three papers presenting a new calculation of the mass 
of the Galaxy based on radial velocities and distances measured for a sample of some 
100 faint 16 < B < 20 field blue horizontal-branch (BHB) stars. This study aims to 
' reduce the uncertainty in the measured mass of the Galaxy by increasing the number 

t^J- , of halo objects at Galactocentric distances r > 30 kpc with measured radial velocities 

[-•««. ■ by a factor five. Faint A-type stars in the Galactic halo have been identified from 

f^*) ' UBjR photometry in six UK Schmidt fields. These samples include field BHB stars as 

CN| . well as less luminous stars of main-sequence surface gravity, which are predominantly 

field blue stragglers. We obtain accurate CCD photometry and spectra to classify 
^ i these stars. This paper describes our methods for separating out clean samples of 

O BHB stars in a way that is efficient in terms of telescope time required. We use the 

I ' high signal-to-noise ratio (S/N) spectra of A-type stars of Kinman, Suntzeff & Kraft 

(published in 1994), and their definitive spectrophotometric A classifications, to assess 
the reliability of two methods, and to quantify the S/N requirements. First we revisit, 
refine and extend the hydrogen line width versus colour relation as a classifier (here 
called the Dq ^-Colour method). The second method is new and compares the shapes 
of the Balmer lines. With this method (here called the Scale width-Shape method) 
there is no need for colours or spectrophotometry. Using the equivalent width of the 
Ca II K line as an additional filter we find we can reproduce Kinman, Suntzeff & 
Kraft's A classifications with both methods. In a sample of stars with strong Balmer 
lines, EW H7 > 13A (equivalent to the colour range < (B - V) Q < 0.2), halo BHB 
stars can be separated from halo blue stragglers reliably. For the spectroscopy (i.e. 

both classification methods) the minimum required continuum S/N is 15 A . For the 
Do.i 5 -Colour method (B — V)o colours accurate to 0.03 mag. arc needed. 
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1 INTRODUCTION 

The mass of the Galaxy is a key quantity for our under- 
standing of the nature and distribution of dark matter. Al- 
though we know a great deal about the contents and prop- 



* Based on observations obtained at the Jacobus Kapteyn Tele- 
scope, the Isaac Newton Telescope, and the William Hcrschel 
Telescope, La Palma, and the Anglo-Australian Telescope, Siding 
Spring Observatory, Australia. 



erties of the halo of the Galaxy our current estimate of the 
total mass is subject to a large uncertainty. In part this 
is due to the Sun's location within the disk of the Galaxy 
which means it is difficult to determine the rotation curve 
accurately outside the solar radius. Beyond ~ 20 kpc mea- 
sures of the enclosed mass of the halo have relied principally 
on the kinematics of satellite galaxies and globular clusters. 
The current state of the art is the analysis by Wilkinson 
& Evans (1999, hereafter, WE99) who calculate the mass 
within 50 kpc to be 5.4+3 g x 1O 11 M , and the total mass 
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to be 1.9li'y x 10 12 Mq. The quoted uncertainties are larger 
than in earlier studies (e.g. Little & Tremaine, 1987, Zaritsky 
et al. 1989, Kochanek, 1996) and were determined by Monte 
Carlo simulation of artificial data sets. The two principal 
sources of error are i) the large uncertainties in the proper 
motions for the six satellites for which measures exist, and 
ii) the small size of the dataset - there are currently only 
27 known satellites at Galactocentric distances > 20 kpc. 
Therefore, we are in the peculiar situation where the mass 
profile of the Galaxy is less well determined than for some 
nearby spiral galaxies. There is no possibility of substan- 
tially increasing the number of known satellite galaxies and 
globular clusters and the best prospect for improving the 
measurement of the mass is through isolating large numbers 
of another distant halo-tracer. Field blue horizontal branch 
(BHB) stars should provide just such a tracer and WE99 
calculate that to reduce the uncertainty on the total mass 
to 20% requires a sample of 200 distant BHB stars. 

This is the first in a series of three papers presenting a 
new calculation of the mass of the Galaxy using radial veloc- 
ities of BHB stars. Field BHB stars are luminous standard 
candles that are abundant in the Galactic halo (e.g. Yanny et 
al. 2000) , and for nearly twenty years, since the study of Pier 
(1983), have presented an under-exploited resource with 
which to measure the density profile and phase space struc- 
ture of the Galaxy halo out to large distances, ~ 100 kpc. 
A number of dynamical analyses of rather small samples of 
BHB stars have been published (e.g. Sommer-Larsen, Chris- 
tensen & Carter, 1989, Norris & Hawkins, 1991, Arnold & 
Gilmore, 1992). Unfortunately, samples of field A-type stars 
in the halo include not only BHB stars but also stars of main 
sequence surface gravity, field blue stragglers, that are some 
2 magnitudes less luminous. Progress towards the goal of 
acquiring a large sample of distant BHB stars has been slow 
because of the difficulty of separating out the BHB stars 
without the investment of large amounts of telescope time. 
In this first paper we describe our procedures for classifying 
samples of halo A-type stars, and present a new efficient 
method that requires only spectroscopic observations of in- 
termediate signal-to-noise ratio. In Paper II we will present 
photometry and spectroscopy of faint 16 < B < 20 candi- 
date BHB stars in two northern high Galactic latitude fields 
and four southern fields. Paper III will contain the dynami- 
cal analysis of the new sample of confirmed BHB stars and 
a new estimate of the mass of the Galaxy. 

The structure of this paper is as follows. In §2 we pro- 
vide a brief outline of the basic parameters of our new survey 
for distant field BHB stars in the halo. §3 contains a sum- 
mary of previous methods employed to sieve out BHB stars 
from samples of halo A-type stars. In §4 we present the de- 
tails of the two classification methods we employ, called the 
Do.is-Colour method, and the Scale width-Shape method. 
We describe the procedures used to measure the Balmer line 
profiles and the equivalent width of the Ca II K line, and 
quantify the random and systematic errors. In §5 we use the 
high signal-to-noise ratio (S/N) spectra of Kinman, Suntzeff 
& Kraft (1994; hereafter KSK) and their spectrophotometric 
classifications to quantify the S/N requirements for apply- 
ing these methods to samples of faint halo stars. In §6 we 
discuss the use of the Ca II K line as a metallicity indicator. 
Finally, §7 provides a summary of the main conclusions of 
the paper. 



2 THE SURVEY 

In a UBV two-colour plot (or its equivalent such as ugr or 
UBjR) of a high Galactic latitude field, halo A-type stars 
are identifiable as a faint continuation of the spectral se- 
quence to hotter types, beyond the main sequence turnoff 
at spectral type F, with colours 0.0 < (B — V)o < 0.2 
(see e.g. Yanny et al. 2000; Fig. 1). As illustrated in strik- 
ing fashion by Yanny et al. using Sloan Digital Sky Survey 
(SDSS) data, these A-type stars include not only luminous 
field BHB stars, of absolute magnitude My ~ 0.7, but also 
stars of main-sequence surface gravity that are some two 
magnitudes less luminous. The nature of these stars is not 
entirely clear but recent work supports the notion that the 
majority are blue stragglers in binary systems created by 
mass transfer when the companion star overfills its Roche 
lobe (Preston & Sneden 2000). As noted by Preston, Beers 
& Shectman (1994) this means that in the field blue strag- 
glers are much more common than BHB stars, whereas in 
globular clusters they are usually rarer. This difference could 
be explained by the destruction of wide binaries in globular 
clusters. This picture is supported by Carrera et al. (2002) 
who find an intermediate value of the ratio of the numbers of 
main-sequence-gravity A stars and BHB stars for the Ursa 
Minor dwarf galaxy, which has a stellar density intermedi- 
ate between that of the field and that of a globular cluster. 
For these reasons in the remainder of this paper we refer to 
these distant halo A-type stars with main-sequence surface 
gravities as field blue stragglers. 

The combination of their substantial luminosity, en- 
suring they can be detected to large distances, and their 
very small spread in absolute magnitude, making them ef- 
fective standard candles, means BHB stars are ideal dynam- 
ical tracers. For the same apparent magnitude the field blue 
stragglers are a factor 2.5 less distant and it is therefore cru- 
cial to separate the two populations. Our survey employs 
APM scans of UK Schmidt Telescope photographic plates. 
We use pairs of plates in each of the U, Bj, and R bands, 
in six fields. The coordinates of the six survey fields are 
provided in Table 1. The fields were selected to give good 
coverage in opposing directions above and below the Galac- 
tic plane, with a range of Galactic longitudes, subject to the 
availability of suitable plate material in the UK Schmidt 
Telescope archive - which effectively restricts the search to 
negative declinations. There was a preference for high Galac- 
tic latitudes to minimize extinction. The different lines of 
sight through the Galaxy provide complementary informa- 
tion on the anisotropy of the stellar orbits. Further details 
of the plate material, observing programme and data pro- 
cessing are provided in Paper II. With plate material from 
the UK Schmidt Telescope A-type stars can be reliably sep- 
arated from F stars down to B ~ 20. We obtain accurate 
CCD BV photometry and spectra at ~ 3A FWHM resolu- 
tion to classify these stars. As described in Paper II we se- 
lect candidate BHB stars by colour in the magnitude range 
16 < B < 20, which corresponds (assuming My ~ 0.7) to 
the distance range 11 — 72 kpc for BHB stars, but to only 
5 — 29 kpc for the field blue stragglers. Pulsating RR Lyrae 
stars have also provided a fruitful source of HB standard 
candles (Hawkins 1984, Ivezic et al. 2000). These stars re- 
side in the instability strip, redward of B-V = 0.2, with AB 
?s 1 mag. and with periods around 0.5 days. However, a num- 
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Field name 


1 


b 


RA 


Dec 








J 2000 


SGP 


250 


-89 


55 


-27 47 


SA94 


175 


-50 


2 53 


12 


F358 


236 


-54 


3 38 


-34 50 


F854 


244 


45 


10 23 


-0 15 


F789 


299 


58 


12 43 


-5 16 


MT 


37 


-51 


22 06 


-18 39 



Table 1. Galactic and equatorial coordinates of the six survey 
fields, ordered by Right Ascension. 

ber of repeat observations are required to isolate these stars, 
making it time consuming to gather observational data. Red 
horizontal branch stars, which are also standard candles, 
have colours (B — V)o > 0.4, similar to the colours of F 
stars, and therefore cannot be identified effectively using 
broadband photometry. 



3 PREVIOUS BHB STAR CLASSIFICATIONS 

To classify A stars the problem is to separate out the gi- 
ants, i.e. the BHB stars, from the dwarfs, i.e. main sequence 
A stars or blue stragglers, using an indicator that depends 
on surface gravity. In samples of relatively nearby stars the 
dwarfs (i.e. main sequence A stars in the disk in this case) 
greatly outnumber the BHB stars, and also the two popu- 
lations have very different kinematic properties. In such cir- 
cumstances, especially for a dynamical study, an extremely 
reliable technique is required to ensure contamination of the 
BHB-sample is minimized, and high-resolution spectroscopy 
is necessary (Kinman et al. 2000). Fortunately, in the outer 
halo, far from the disk, the number of BHB stars and po- 
tential contaminants (here blue stragglers) are comparable. 
Furthermore the two populations will have much more simi- 
lar kinematic properties, so the problem is somewhat less de- 
manding. Our goal is to identify an efficient method for the 
classification of BHB stars in the halo that is largely com- 
plete while also ensuring the contamination from non-BHB 
stars is 10%. Under these requirements high resolution 
spectroscopy is not necessary. 

The study by KSK of the problem of classifying distant 
A-type stars in the Galactic halo has provided the bench- 
mark for future work. KSK studied three surface gravity 
indicators, each as a function of (B — V)o colour: i) the 
strength of the Balmer jump, which may be quantified by 
the Stromgren index ci, or some similar close equivalent, ii) 
the steepness of the Balmer jump, quantified by the spec- 
trophotometric A parameter, iii) the width of the Balmer 
lines. (Note that broadband C/i3V-photometry on its own 
provides only a crude separation of the populations and was 
not considered.) Of the three methods the A parameter pro- 
vides the cleanest separation and is the preferred indica- 
tor. However, measurement of the A parameter is impracti- 
cal for faint stars. A S/N~ 40 A ^ is required, which, at 
V = 19, takes some two hours of integration on a 4m tele- 
scope in dark time. Furthermore, a wide slit is required for 

' In this paper all values of S/N are quoted as A by dividing 
the S/N per pixel by the square root of the pixel size in A. 



spectrophotometry, necessitating the acquisition of a second 
spectrum, with a narrower slit, to measure the star's radial 
velocity accurately. Method i) is also impractical for very 
similar reasons. 

The third method, based on the width of the Balmer 
lines as a function of (B — V)o colour, is the most efficient, 
since the spectrum also provides the radial velocity. This 
method has been used in all previous studies of the dynam- 
ics of faint BHB stars. Flynn, Sommer-Larsen & Christensen 
(1994) classified a sample of A-type stars using the Strom- 
gren index c\ and demonstrated that the stars classified as 
BHB also cleanly separate out in a plot of line-width against 
(B — V)o colour. Similarly, KSK found that stars classified 
BHB by the A parameter separate from blue stragglers on 
the basis of line-width. However, neither Flynn et al. nor 
KSK quantified the spectroscopic S/N required or the frac- 
tional contamination of BHB samples at a given S/N. Fi- 
nally, Wilhelm, Beers & Gray (1999a) made a detailed study 
of the problem, adding U photometry to the BV photome- 
try and line-widths. Nevertheless, although the U — B colour 
is a useful indicator of surface gravity for redder colours 
0.2 < (B — V)o < 0.4, it adds less information over the 
range 0.0 < (B - V) < 0.2. 



4 NEW TECHNIQUES FOR CLASSIFYING 
FAINT HALO A-TYPE STARS 

In this section we describe new techniques for classifying 
faint halo A-type stars efficiently. In considering the best 
way to separate BHB stars and blue stragglers we first re- 
visited the use of the line width versus colour relation as a 
discriminant. We have made a number of small refinements 
over the work of KSK: i) we use a functional fit to mea- 
sure the line widths, removing subjectivity and allowing us 
to quantify the uncertainties, ii) we use a slightly different 
measure of line width than KSK, called D0.15, iii) we use 
the improved reddening corrections of Schlegel, Finkbeiner 
& Davis (1998), and iv) we use the Ca II K absorption line 
EW as an additional classification filter. In the rest of this 
paper we refer to this method as the Do.is-Colour method. 

In the course of this work we discovered that the depen- 
dence of the profile (i.e. the detailed shape) of the Balmer 
lines on both temperature and surface gravity is measurable. 
This led to what we have called the Scale width-Shape classi- 
fication method. With this method there is no need to obtain 
colours or spectrophotometry. These two classification meth- 
ods are the subject of the remainder of the paper. Note that 
the two methods are not entirely independent since they 
both use parameters derived from a fit to the line shape. In 
§5 we consider the circumstances under which one or other 
is more useful. 

In this section we summarise details of the KSK-sample 
of spectra used in developing the classification methods. We 
then describe the procedure for measuring the profiles of 
the Balmer lines, and its application to the two methods of 
classification. Lastly we explain the procedure for measuring 
the Ca II K line EW, and the uses of this measure. 



© 0000 RAS, MNRAS 000, 000-000 



4 L. Clewley, et al. 



4.1 The KSK data 

The high S/N spectra of stars analysed by KSK were kindly 
made available in electronic form by Dr Nick Suntzeff. The 
total sample includes 214 stars ranging in spectral type from 
B to F. The majority of the sample are A-type stars, includ- 
ing BHB stars, main-sequence A stars, and blue stragglers. 

The methods of classification we apply are effective in 
the temperature range where the Balmer lines are strong and 
do not work blueward of (B — V)o = 0. Because the Scale 
width-Shape method does not use colours, but both methods 
use spectroscopy, we need a spectroscopic criterion, rather 
than a colour criterion, for defining the type of star to which 
the classification methods apply. We select all stars with 
EW H7 > 13 A. This sample contains 131 stars. We refer to 
this as the 'KSK total sample'. The EW limit corresponds 
approximately to the colour range < (B — V)o < 0.2. The 
EW limit was chosen by trial and error to produce the largest 
clean samples of BHB stars. With a lower limiting EW we 
include stars with (B - V)o < and (B - V)o > 0.2 where 
discrimination is more difficult. With a higher limiting EW 
we remove BHB stars from the sample but retain most of 
the blue stragglers. 

The main-sequence stars in this sample are typically of 
solar metallicity. Since we are interested in classifying lower- 
metallicity halo stars we focus in particular on a subsample 
of fainter high-Galactic latitude stars, consisting of 66 stars, 
of magnitudes 13.0 < V < 16.5 in two fields, SA57 at the 
Northern Polar Cap, and RR7 in the direction of the Galac- 
tic anticentre. These stars have S/N in the range 23 to 80 
A , with a mean of 66 A . We refer to this subsample as 
the 'KSK halo sample'. The KSK total sample is useful for 
showing the broad trends in the classification parameters, 
because with the larger sample the sequences of the two 
populations are better defined (e.g. Fig. but because of 
the issue of metallicity we use the KSK halo sample only in 
defining the classification boundaries in plots of the relevant 
parameters. 

In comparing line widths measured for the same star 
in the KSK spectra observed in different runs, we discov- 
ered a small systematic difference between the spectra in 
the two halves of the dataset. The origin of the difference 
was traced, with the assistance of Dr Suntzeff, to the details 
of the flux-calibration procedures employed for different ob- 
serving runs. We were able to correct the half of the data in 
error by comparing the KSK spectra of the flux standard, 
flux calibrated by itself, to the original Massey et al. (1988) 
spectrum. 



4.2 Functional fit to the Balmer lines 

The results of a functional fit to the Balmer lines are used in 
both classification methods. The width Do. 15 is determined 
from the fit, while the parameters of the fit themselves are 
used in the Scale width-Shape method. 



4-2.1 Continuum fit 

Before fitting the Balmer lines each spectrum is normalised 
to the continuum by fitting a polynomial of degree three ^] to 
regions of continuum well away from the wings of the lines. 
Because the Balmer lines are so broad care must be taken in 
fitting the continuum, and we spent some time experiment- 
ing with the degree of the polynomial and with different 
wavelength ranges for fitting. The intervals chosen for the 
fit were 3863-3868, 3902-3925, 4020-4048, 4146-4275, 4388- 
4494 A. While the final choice of procedure is inevitably sub- 
jective this is not a concern provided other observers who 
follow the same procedure obtain the same results, within 
the errors. This will be the case provided there is no system- 
atic trend of the measured parameters with S/N. 

To test for systematics we created artificial spectra by 
adding noise to high quality spectra of four stars. For each 
star we created 1000 spectra of a specified S/N, for sev- 
eral values of S/N in the range 7 to 30. At a given S/N 
we measured the lines in each artificial spectrum, and then 
calculated the mean and the error on the mean for each 
measured parameter. For no parameter, at any S/N, was 
the mean value of the parameter inconsistent with the value 
measured for the original high S /N spectrum. In other words 
we found no significant systematic errors associated with the 
continuum fit. 

4.2.2 Sersic profile fit 

To fit the line profiles in the normalized spectra we ex- 
perimented with several functions including an exponen- 
tial, a —5/2 power law (Pier, 1983), a general power-law, 
a Lorentzian, a Gaussian (Beers, et al. 1992), and a Voigt 
profile (Wilhelm et al. 1999b). Finally we adopted the Ser- 
sic profile (Sersic, 1968), which is the exponential of a power 
law: 

(X — Xq \\ c ~ 
6 J 

This is the function with the fewest free parameters that 
gives a good fit while remaining well behaved. Because the 
lines have a sharp core, in fitting we convolve model spectra 
with a Gaussian of FWHM equal to the instrument resolu- 
tion i.e. we measure the deconvolved profile. The function 
has four free parameters. The parameter xo is the wave- 
length of the line centre. The parameter a is the depth at 
the line centre. The parameter b is the scale length of the 
exponential, and so gives a measure of the line width. We 
refer to b as the scale width in this paper, to distinguish it 
from the other line width measure we use, Do.is- Finally the 
parameter c quantifies the line shape. Note that c — 2 corre- 
sponds to a Gaussian, c = 1 to an exponential, and c — 1/4 
to the de Vaucouleurs profile (de Vaucouleurs 1948). For 
the same value of b, a smaller value of c produces broader 
wings to the line. An example fit for a typical BHB star in 
the KSK sample is shown in Fig. [j] (b). The shapes of the 
Balmer lines are discussed further in §4.2.4. 

We found that the parameter a is not useful for sepa- 
rating the two populations, but that in a plot of b against 
c the two populations are distinct (Fig. Q). Because of 

t in IRAF parlance this is a polynomial of order four 
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Figure 1. (a) Illustration of the difference in line profile between 
BHB stars and blue stragglers. The solid line plots the intrinsic 
(i.e. deconvolved) profile of the star in (b) with b = 8.9A and 
c = 0.91. The dashed line shows a blue straggler with the same 
temperature c = 0.91 and a broader scale width, b = 12A. Al- 
ternatively, the dotted line shows a blue straggler which has the 
same scale width as the BHB star, b = 8.9A, but is cooler, c = 0.7, 
and has broader line wings, (b) Convolved fit of the H5 Balmer 
line for a typical BHB star in the KSK dataset. This star has 
measured (deconvolved) parameters b = 8.9A, c = 0.91. 



the relatively low resolution of the KSK spectra (typical 
FWHM= 4 A), and the sharpness of the line core, the (de- 
convolved) parameter a is difficult to measure accurately. 
This is unfortunate because errors on the interesting param- 
eters b and c are correlated with errors in a. The parameter 
a is approximately constant for the strong-lined stars in the 
KSK total sample. For this reason we decided to fix a at 
the average value a = 0.83. This is satisfactory for the KSK 
data, and our own programme spectra, but might not pro- 



Line 


Rcstframe wavelength (A) 


Bandpass (A) 


Ca II K 


3934 


3919-3949 


H<5 


4102 


4032-4172 


Ho 


4341 


4271-4411 



Table 2. Wavelength regions used in fitting the profiles of the 
Balmer lines and measuring the Ca II K EW. 



vide an adequate fit for higher resolution data with higher 
S/N. 

To summarize, in measuring the line profiles we com- 
pute the minimum-x 2 fit of the following function convolved 
with a Gaussian of FWHM equal to the spectral resolution: 



- x 



y = 1.0 - 0.83 exp 



We found satisfactory fits to all the stars in the KSK 
total sample with this profile. (This modified function, with 
a fixed at 0.83, is not appropriate for B or F stars however.) 
The Scale width-Shape method uses the parameters 6 and c 
determined in this way. It should be noted that the separa- 
tion boundary would be slightly different for measurements 
where a is a free parameter 

The measurements were restricted to the and HS 
lines. The He line is blended with the Ca II H absorp- 
tion line, while the higher-order Balmer lines crowd to- 
gether and it becomes difficult to define the continuum. 
The KSK observations do not cover the H/3 line. The wave- 
length intervals used in making the fit to the H'y and 
HS lines are provided in Table 2. In making the fit dis- 
crepant pixels are iteratively a— clipped. The measured line 
width may depend on metallicity due to the presence of 
weak metal lines that are not clipped out. In an attempt 
to reduce the effect of metallicity we have identified the 
10 strongest metal lines that lie within the profile of ei- 
ther the Hj or HS line. These are listed in Table |^. Pix- 
els at the doppler-shifted wavelength of the metal lines are 
masked out. We return to the issue of line blanketing below. 
The fitting routine runs within IRAF and is available at 
http : //astro . ic . ac . uk/Research/extragal/milkyway . html 



4-2.3 Do. 15 -Colour method 

Over the colour range < (B — V)o < 0.2 the Balmer 
lines are strong and broad, and A stars can be separated 
according to surface gravity by plotting line width (suitably 
defined) against colour, as shown for example in Fig. ^. 

In order to measure the Balmer line-widths most au- 
thors have followed Rodgers, Harding & Sadler (1981) and 
used Do. 2, defined as the width measured at a flux level of 
0.8 of the continuum flux level. KSK used a similar measure 
but the line-depth was determined relative to the bottom of 
the Balmer line rather than from the level of zero flux. This 



3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation. 
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Line 


Wavelength (A) 


Fe I 


4045.8 


Sr II 


4077.7 


Fe I 


4173.0 


Ca I 


4226.7 


Fe I 


4271.8 


Ti II 


4300.1 


Fe II 


4351.8 


Fe II 


4385.4 


Fe II 


4416.8 


Mg II 


4481.3 



Table 3. Restframe wavelengths of metal lines masked out by 
the fitting routine. 

definition produces a measure of the width closer to the con- 
tinuum, fully exploiting the gravity difference in the wings. 
The advantage of this is made clear by the curves plotted in 
Fig. |l| (a). However because the core of the Balmer lines is 
so sharp the (convolved) line depth depends on the spectral 
resolution. To combine the best features of the two different 
measures we have defined Do. 15 which is the width measured 
at a flux level 0.85 of the continuum flux level. 

In order to measure the parameter Do. 15 using the Ser- 
sic profile fit, the profile is re-expressed in terms of the three 
free parameters Do. 15, %o, and c. (For reference, the relation 
between the parameters is D .i5 = 2&ln(0.83/0.15) 1/c .) The 
la uncertainty on Do. 15 , marginalising over the other param- 
eters, was computed by determining the increment AD0.15 
that produced A^ 2 = 1.0 when minimizing on the other two 
parameters. We checked the errors by two methods. First we 
measured the scatter in repeat measurements of the same 
star. Second we took high S/N spectra and created random 
realizations of lower S/N spectra, and measured the scat- 
ter in the fitted parameters. In both cases we found good 
agreement with the errors output by the fitting routine. 

4-2.4 Scale width-Shape method 

By plotting different combinations of two of the four param- 
eters a, b, c, and (B — V)o for the KSK sample we found 
that the parameters b and c provide an alternative method 
for distinguishing between BHB stars and blue stragglers. As 
shown in Fig. ^ over the colour range 0.0 < (B — V)o < 0.2 
the parameter c, which quantifies the shape of the line, is 
tightly inverse-correlated with colour, following the linear 
relation c = 0.93 — 0.98(D — V)o. Evidently c provides a 
measure of temperature for A-type stars. At a fixed value 
of c (temperature) stars of higher surface gravity have larger 
values of the scale width b. In a plot of b against c (e.g. Fig. 
BHB stars separate from blue stragglers. 

The differences between the line profiles for the two dif- 
ferent types of star are illustrated in Fig. |l| (a). Blue strag- 
glers have larger scale widths than BHB stars of the same 
temperature. Alternatively blue stragglers of the same scale 
width as a BHB star are cooler (smaller c), and therefore 
have sharper line profile cores and broader line wings. 

The Scale width-Shape method is effective for stars with 
strong Balmer lines i.e. A-type stars. A major advantage of 
this method of discrimination lies in the fact that no pho- 
tometric measurements are required to separate the popu- 
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Figure 2. Plot of the shape parameter c (averaged for the H<5 and 
H7 lines) against (B — V)o for the KSK halo sample, i.e. A— type 
stars with Balmer EW> 13 A. Evidently c provides a measure of 
temperature. 

lations. As b and c are correlated parameters we compute 
the error ellipse, marginalizing over the other parameter, 
xO. In plotting the uncertainties we show semi-major and 
semi-minor axes of the ellipse, marking the 68% confidence 
interval for each axis in isolation, i.e. marginalizing along 
the other axis. We checked the errors in the same way as for 
the D . 15 -Colour method and again found good agreement. 

4.3 Combining US and H7 line widths 

For each parameter we have measurements for the two lines 
H<5 and H7 in each spectrum. These can be combined to re- 
duce the errors. For each of the width parameters Do. 15 and 
b we found small but significant differences between the val- 
ues measured for the H<5 and the H7 lines, with a dependence 
of the difference on temperature. This effect is illustrated in 
Figs ^| and ||. The first of these plots the difference in Do. 15 
between the two lines against (B — V)rj. The second of these 
plots the difference in b between the two lines against the 
parameter c. We have chosen the H7 line as our reference, 
and we have used the linear fits to calibrate the values of 
Do. 15 and b from H<5 to H7 and then calculated weighted 
averages. For the parameter c we found no such trend, and 
so we simply computed the weighted average of the c-values 
for the two lines. 

We have not discovered a convincing explanation for the 
trends seen. We suspect it is an effect of line blanketing from 
weak metal lines, which would explain the temperature de- 
pendence. However if this is the case one would expect there 
to be a correlation between the residual from the fit and the 
metallicity, as estimated from the Ca II K line. However we 
found no significant correlation in a plot of these two quan- 
tities. 
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Figure 3. Difference between Do. 15 for the H7 and Hc5 lines 
as a function of colour, for the KSK halo sample. A systematic 
difference in line width is evident. A linear least-squares fit to the 
data (with la errors) gives Do. 15 (7) — A). 15 (S) = (0.44 ±0.16) — 
(6.42 ± 1.17)(B - V) , shown as the solid line. The dotted line 
Do. 15(7) — Dq. 15(5) = 0.0 is also shown as a guide. 




0.6 0.7 0.8 0.9 1 

(c( 7 <5)) 

Figure 4. Difference between 6(7) and b(8) as a function of c, 
for the KSK halo sample. A systematic difference in line width 
is evident. A linear least— squares fit to the data (with lcr errors) 
gives 6(7) - b(S) = (1.71 ± 0.33) - (1.31 ± 0.40)(c( 7 <5)) 



4.4 Ca II K-line EW measurements 

The Ca II K line is the strongest measurable metal line over 
the wavelength range covered by the KSK spectra, and the 
only useful indicator of metallicity in spectra of low S/N, 
such as for the stars in our radial velocity programme. With 
some exceptions (including the A metallic (Am) and pecu- 



liar (Ap) stars) the strength of the Ca II K line at constant 
temperature can be used as a reliable indicator of the metal- 
licity of A-type stars (e.g. Pier 1983, Beers et al. 1992, KSK). 
We measure the EW of the Ca II K line by minimum— \ 2 
fitting a Gaussian to the continuum divided spectrum over 
the wavelength range given in Table 2. This line is much 
weaker than the Balmer lines, and to reduce the error the 
central wavelength is fixed at the redshift determined from 
the Balmer lines. 

In §6 we use a plot of EWc a versus (B — V)o to estimate 
the metallicities of the KSK stars. This has two purposes. 
First it allows anomalous high-metallicity halo stars to be 
identified. Secondly knowledge of the metallicity can be used 
to reduce the uncertainty in the absolute magnitudes of the 
BHB stars, and consequently, the errors in the distances (e.g. 
Clementini et al. 1995). The issue of measuring distances is 
deferred to Paper II. 



5 RESULTS: THE BALMER LINES 

In this section we show how the two classification methods 
work in practice. Using measurements of the Balmer lines 
for the KSK sample we present plots illustrating the two 
methods. We then quantify the S/N required to apply the 
classification methods to samples of faint halo A-type stars. 

5.1 Do. 15 -Colour method 

The Dq. 15-Colour method is illustrated in Fig. |B| which plots 
Do. 15, averaged for the two lines, against (B — V)o colour. 
We have used KSK's classification using the spectrophoto- 
metric A parameter to separate the samples into two pop- 
ulations. Note that because this method uses the shape of 
the continuum it is independent of classification that uses 
the line widths. In the two plots stars classified BHB are 
plotted as filled symbols and other stars i.e. blue stragglers 
or main-sequence A stars, hereafter 'A/BS', are shown as 
open symbols. The left-hand plot is for the KSK total sam- 
ple, which includes stars with a range of metallicities, while 
the right-hand plot shows the KSK halo sample only. 

The two populations BHB and A/BS are cleanly sep- 
arated in these plots. The two sequences stand out more 
clearly in the left-hand plot because of the larger number 
of objects. The curves overplotted are derived from model 
spectra. The calculation of the models is described in the Ap- 
pendix. Curves for three surface gravities log g = 3.0, 3.5, 4.0 
are shown, the line width increasing with surface gravity. 
For each value of log g curves for two different metallicities 
[Fe/H] = -2.0, -1.0 are plotted. The right-hand curve of each 
pair corresponds to the higher metallicity. While the curves 
follow the broad trends in the data, in detail they provide 
a rather poor fit. Main sequence stars in this temperature 
range have logg ~ 4.5 whereas the curves with logg = 4.0 
provide the best fit to A/BS sequence near (B — V)o = 0.1, 
while falling below the data at redder colours. 

Since we are interested in classifying samples of faint 
halo A-type stars we concentrate on the right-hand plot in 
Fig. ^ We have chosen the curve computed for logg = 3.5 
and [Fe/H] = -1.0, marked bold, as the classification bound- 
ary. Above the line there are 31 stars classified A/BS by 
KSK and 2 stars classified BHB. Below the line there are 
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Figure 5. Separation of the KSK sample using the Do.i5-Colour method: (a) the KSK total sample (131 stars), (b) the KSK halo 
sample (66 stars). BHB stars, classified as such by KSK using their A method, are marked by filled circles, and stars classified A/BS are 
marked by open circles. The curves show results for model spectra from Kurucz (1993) for log g = 3.0, 3.5, 4.0. Each gravity has two 
lines of metallicity [Fe/H] = -2 (left) and -1 (right). The line log g = 3.5, [Fe/H] = -1, plotted bold, is our chosen classification boundary. 



32 stars classified BHB and only 1 star classified A/BS. For 
reference, in the left hand plot above the line there are 77 
stars classified A/BS by KSK and 3 stars classified BHB, and 
below the line there are 47 stars classified BHB and 4 stars 
classified A/BS. Clearly with data of this quality highly reli- 
able classification is possible with the Do. 15 -Colour method. 
The nature of the few stars below the line classified A/BS 
is discussed in §6. Most appear to be anomalous by virtue 
of their high metallicity and can be identified and removed. 



5.1.1 S/N requirements 

We have used Monte-Carlo methods to assess the spec- 
troscopic and photometric S/N requirements for classifying 
samples of more distant halo A-type stars using the Do. 15- 
Colour method. We take a figure of ~ 10% as the practical 
acceptable limit of contamination of any sample of halo BHB 
stars by blue stragglers. We started with the simplifying as- 
sumption that the KSK halo sample is representative of the 
halo population at fainter apparent magnitudes in terms of 
the ratio of blue stragglers to BHB stars. Then, adopting 
specific values of spectroscopic and photometric S/N, for 
each star in the KSK halo sample we drew a value of Do. 15 
and (B — V)o from the error distributions to create a syn- 
thetic halo sample. By creating thousands of samples and 
counting the number of stars originally classified A/BS or 
BHB that appear below the classification boundary in each 
sample we have estimated the contamination and complete- 
ness of samples of faint halo BHB stars as a function of S/N. 

Fig. |^ shows the results of these simulations. The up- 
per plot shows the percentage completeness of BHB samples 
as a function of spectroscopic S/N, for two fixed values of 
photometric S/N. The lower plot shows the percentage con- 



tamination of BHB samples. For our faint programme stars 
16 < B < 20 we have found it difficult to achieve errors 
in (B — V)o smaller than 3%. For this level of photometric 
accuracy, Fig. ^ shows that a spectroscopic S/N of 15 A -1 
would achieve ~ 87% completeness with contamination by 
blue stragglers at a level of ~ 7%. We consider this accept- 
able. 

Coincidentally this S/N requirement matches the re- 
quirement for measurement of the radial velocity. In our ex- 
perience with spectra of S/N=15A _1 the accuracy of the ra- 
dial velocity measurement is 15 km s~\ which is satisfactory 
considering that the one-dimensional velocity dispersion in 
the halo is ~ 100 km s _1 . 

To summarize, with (B — V)o colours accurate to 3%, 
with spectra of S/N=15 A" 1 , samples of halo BHB stars that 
are ~ 87% complete, and with only ~ 7% contamination 
by blue stragglers, can be compiled using the Do.is-Colour 
method. 

It would be possible to reduce the contamination at the 
expense of sample size, by restricting the sample to a nar- 
rower range of (B — V)o where the discrimination is cleaner. 
For example if we now limit the analysis of the KSK halo 
sample to stars in the range 0.05 < (B — V)o < 0.15 the 
number of BHB stars is reduced from 32 to 19. For this sam- 
ple, with (B — V)o colours accurate to 3% and with spectra 
of S/N=15A _1 , compared to the sample with the broader 
colour range the completeness is increased from ~ 87% to 
93% and the contamination is reduced from 7% to 5%. 



5.2 Scale width-Shape method 

The Scale width-Shape method is illustrated in Fig. |^, which 
plots b against c, averaged for the two lines. As before, the 
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virtue of their high metallicity and can be identified and 
removed. 



5.2.1 S/N requirements 

We have followed the same Monte-Carlo procedure used for 
the Do.is-Colour method to establish the S/N requirements 
for the Scale width-Shape method. The results are plotted 
in Fig. |§| Here the only variable is spectroscopic S/N. The 
completeness and contamination for this method are slightly 
worse than for the Do. 15 -Colour method. For a spectroscopic 
S/N of 15 A -1 the completeness is ~ 82% (compared to 
87%) and the contamination is ~ 12% (compared to 7%). 

We see that with spectroscopy alone we can separate 
the two populations almost as well as with the Do.i5~Colour 
method. Because accurate photometry is not required there 
will clearly be circumstances where the Scale width-Shape 
method is preferred. For example if the photometric accu- 
racy of the parent catalogue from which the A-type stars are 
selected is significantly worse than 3% the colours are not 
useful for classification. It may then be most efficient to sim- 
ply obtain spectra of the candidate BHB stars. This saves 
the time required to obtain accurate photometry, with the 
only penalties a small reduction in completeness and a small 
increase in contamination. This may well be the best strat- 
egy for exploring the outer reaches of the Galaxy's halo using 
the SDSS dataset; stars at 100 kpc distance have g' ~ 21, 
where the error on the Sloan g' — r' colour is 5% (Stoughton 
et al., 2002). 



Figure 6. Monte Carlo simulations of the halo sample of stars 
shown in Fig. |i|(b): (a) completeness of the BHB sample as a 
function of spectroscopic S/N, for 2% and 3% errors in the colours 
(b) contamination of the sample by blue stragglers. 



6 RESULTS: THE CA II K LINE 

In this section we establish the accuracy with which metal- 
licities can be measured using the Ca II K line EW. We then 
show that the estimated metallicities are useful for identify- 
ing interlopers in the BHB samples. 



left-hand plot is for the KSK total sample, while the right- 
hand plot shows the KSK halo sample only. Again the two 
populations BHB and A/BS are clearly separated in these 
plots. The curve shown is the empirical selection bound- 
ary that we have chosen for this method. Unfortunately the 
model curves are a very poor match to the data and there- 
fore have not been plotted. The problem appears to be with 
the shapes of the lines more than with than the widths. We 
find that the model spectra reproduce the slope of the rela- 
tion between c and (B — V)o reasonably well, but that the 
locus for the model spectra is offset to higher values of the 
parameter c. In other words the Balmer lines of the models 
are not spikey enough. 

In the right-hand plot, showing the halo sample, above 
the line there are 31 stars classified A/BS by KSK and 
classified BHB. Below the line there are 34 stars classified 
BHB and only 1 star classified A/BS. For reference, in the 
left hand plot above the line there are 73 stars classified 
A/BS by KSK and stars classified BHB, and below the 
line there are 51 stars classified BHB and 7 stars classified 
A/BS. Therefore the Scale width-Shape method appears to 
be similarly effective in separating the two populations. The 
nature of the few stars below the line classified A/BS is 
discussed in §6. Again, most appear to be anomalous by 



6.1 Measurement of metallicity 

KSK plotted the Ca II K line EW against (B — V)o for stars 
for which independent accurate metallicities are known, to 
determine empirically curves of constant metallicity in this 
parameter space. Instead we have used the theoretical curves 
of Wilhelm et al. (1999a) measured from synthetic spectra. 
In Fig. we plot their isoabundance contours for [Fe/H] 
= -1, -2 and -3 (the solar- metallicity line is explained be- 
low). Metallicities can be determined from this plot using 
measurements of EWca and (B — V)o, by interpolation. To 
assess the accuracy of this plot we have used it to measure 
metallicities for all the BHB stars in the KSK total sam- 
ple for which independent accurate metallicities are known. 
The stars used, listed in Table 4 and plotted in Fig. are 
the BHB stars in M92 ([Fe/H] = -2.2), in M3 ([Fe/H] = 
-1.5), and a sample of nearby field BHB stars (individually 
labeled with the value of [Fe/H]). The accurate metallicities 
are taken from KSK and are listed in Table 4 in the column 
headed [Fe/B] K . 

Our interpolated measured metallicities are listed in the 
last column of Table 4, labeled [Fe/H]o. There is quite good 
agreement between our interpolated estimates and the accu- 
rate determinations. For the sample of 18 stars we measure 
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Figure 7. Separation of the KSK sample using the Scale width-Shape method: (a) the KSK total sample (131 stars), (b) the KSK halo 
sample (66 stars). BHB stars, classified as such by KSK using their A method, are marked by filled circles, and stars classified A/BS 
are marked by open circles. Curves showing results for model spectra provide a very poor fit and have been omitted. The dashed curve 
indicates the adopted classification boundary. 



a mean difference [Fe/H]c-[Fe/H]x = 0.10 with standard 
deviation ±0.24. For the three sub-samples M92, M3, and 
the field BHB stars the results are 0.21 ± 0.20 (8 stars), 
-0.10±0.17 (3 stars), 0.06±0.27 (7 stars) respectively. The 
scatter of the points could arise from a variety of sources 
including the difficulty in measuring the continuum, con- 
tamination of the line by interstellar absorption, and uncer- 
tainties in the accurate estimates themselves (see KSK for a 
discussion). Therefore in using this plot to measure metallic- 
ities of other stars we add an error of 0.3 dex in quadrature 
to the random error, based on the measured scatter. Addi- 
tional scatter may be introduced by the spectral peculiarities 
of the Am and Ap stars (discussed in KSK, Wilhelm et al. 
1999a, 1999b). 

To measure the metallicities of other A-type stars over 
the full range of metallicities we have defined also a solar- 
metallicity contour by fitting a straight line to the data 
for main-sequence A stars in the Pleiades and Coma clus- 
ters, plotted as open triangles. Because the lines are con- 
verging towards bluer colours we consider this plot reli- 
able for estimating metallicities only for colours redder than 
(B - V)o > 0.05. 



6.2 Identification of interlopers 

We return now to the issue of interlopers i.e. the few halo 
stars that are classified by our methods as BHB but classified 
by KSK as A/BS. We have discovered that most of these 
interlopers have anomalously high metallicity. 

If we consider first the -D0.15 -Colour method, then for 
the KSK halo sample, shown in Fig. ^ (RHS), there is one in- 
terloper amongst the 33 stars below the classification bound- 
ary, the star numbered RR7-70 by KSK. This star is the only 



star below the line with estimated metallicity [Fe/H]> —0.5. 
In the KSK total sample, Fig. ^ (LHS), there are 4 inter- 
lopers amongst the 51 stars below the classification bound- 
ary, and similarly they are the only stars with metallicity 
[Fe/H]> -0.5. 

Turning to the Scale width-Shape method, in the halo 
sample, Fig. ^ (RHS), there is one interloper amongst the 
33 stars below the classification boundary, the star RR7-70, 
the same interloper in Fig. ^| (RHS). This star is the only 
star below the line with estimated metallicity [Fe/H]> —0.5. 
In the total sample, there are 7 interlopers amongst the 58 
stars below the classification boundary. For three of these 
we are unable to estimate metallicities reliably as they have 
colours (B — V)o < 0.05. Three of the remaining four are 
the only stars with metallicity [Fe/H]> —0.5. 

We conclude that there is a very small proportion of 
stars classified BHB by our method that are probably not 
BHB, and that these can be identified by their anomalously 
strong Ca II K line at the given temperature. Applying this 
additional filter with high S/N data, samples of BHB stars 
that are almost perfectly clean and complete can be pro- 
duced using the two classification methods discussed in this 
paper. 



7 SUMMARY 

The purpose of this paper has been to set out efficient meth- 
ods for classifying reliably halo A-type stars into the classes 
BHB and blue straggler, and to quantify the S/N require- 
ments. We disregarded the spectrophotometric A method 
and Stromgren u photometry as possible methods of classifi- 
cation because they require too much telescope time, bearing 
in mind that additional spectroscopy is required to measure 
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Figure 8. Monte Carlo simulations of the halo sample of stars 
shown in Fig. 0(b): (a) completeness of the BHB sample as a 
function of spectroscopic S/N, for 2% and 3% errors in the colours 
(b) contamination of the sample by blue stragglers. 



radial velocities. Instead we have concentrated on measure- 
ments of the Balmer lines H7 and H5 from intermediate- 
resolution (~ 3A) spectroscopy (not spectrophotometry). 
The Do. 15 -Colour method plots the width of the Balmer 
lines against (B — V)o colour. The main improvements we 
have introduced over previous applications of this method 
are profile fitting of the Sersic function, and quantification 
of the errors. We find that with spectra of S/N= 15 A -1 
and colours accurate to 3% samples of halo BHB stars will 
be ~ 87% complete with only ~ 7% contamination by blue 
stragglers. Spectra of this S/N provide radial velocities ac- 
curate to 15 km s _1 . The small contamination of the sample 
can be further reduced by identifying stars with anomalously 
strong Ca II K absorption for their temperature. 

A second classification technique, the Scale width-Shape 
method, compares the shapes of the Balmer lines by plot- 
ting two of the parameters of the profile fit. We find that 
with spectra of S/N= 15 A -1 samples of halo BHB stars 
will be ~ 83% complete with only ~ 12% contamination by 
blue stragglers. In other words without the need for colours 
the samples are almost as clean and complete as with the 
Axis -Colour method. This argues that if accurate colours 
< 3% are available for the sample the best method is the 
Do.i5~Colour method. However where the accuracy of the 
photometry is worse the most efficient method for studying 
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Figure 9. Ca II K line (3933 A) EW for selected stars in KSK, 
listed in Table 4. The curves represent lines of metallicity for 
[Fe/H] = -1.0, -2.0 and -3.0 taken from Wilhelm et al. (1999a). 
The straight line represents a best fit to stars in the Pleiades and 
Coma clusters (open triangles) assumed to be of solar metallicity. 
Filled squares are stars from M92 ([Fe/H] = -2.2) and filled tri- 
angles are from M3 ([Fe/H] = -1.5). Field BHB stars are shown 
by open circles labeled with the KSK adopted values. 
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EWca 




(B- 


- V) 


[Fe/H] K 


[Fe/H] 


M92 11-23 


0.50 





07 





06 


-2 


2 


-2.1 


M92 IV-27 


1.03 





07 





15 


-2 


2 


-1.9 


M92 S-20 


0.65 





07 





13 


-2 


2 


-2.3 


M92 S-24 


0.52 





05 





01 


-2 


2 


-1.7 


M92 VI-10 


0.64 





07 





11 


-2 


2 


-2.2 


M92 XII-01 


1.05 





05 





17 


-2 


2 


-2.0 


M92 XII-09 


0.99 





07 





13 


-2 


2 


-1.8 


M92 XII-10 


0.53 





06 





05 


-2 


2 


-1.9 


M3 182 


1.07 





07 





08 


-1 


5 


-1.5 


M3 11-11 


4.05 





13 





42 


-1 


5 


-1.8 


M3 VI-18 


1.63 





06 





16 


-1 


5 


-1.5 


HD 2857 


1.76 





06 





17 


-1 


7 


-1.5 


HD 14829 


0.41 





05 





01 


-2 


2 


-1.9 


HD 60778 


1.11 





02 





08 


-1 


4 


-1.4 


HD 74721 


0.55 





01 





02 


-1 


3 


-1.7 


HD 86986 


0.92 





02 





09 


-1 


7 


-1.7 


HD 109995 


0.64 





01 





05 


-1 


7 


-1.8 


HD 161817 


1.61 





02 





13 


-1 


8 


-1.4 



Table 4. Ca II K line EW and the calculated [Fe/H] ([Fe/HlrO 
compared to those in KSK ([Fe/H]x) for stars shown in Fig. M. 



the dynamics of the Galactic halo would be the Scale width- 
Shape method, saving the need to obtain accurate photom- 
etry. This may be the best strategy for studying the outer 
halo of the Galaxy beyond r = 100 kpc, with candidate BHB 
stars selected from the SDSS database. 
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APPENDIX A: MODEL SPECTRA 

We assembled a grid of spectral line profiles using LTE 
model atmospheres created by Kurucz (1993) and the spec- 
trum synthesis code SYNSPEC (Hubeny et al. 1994). We 
created spectra over the wavelength range 3500 — 4500 A. 
Detailed Stark broadening tables were used for the hydro- 
gen lines and the flux distribution was calculated every 
0.01 A. We examined models with logg = 3.0,3.5,4.0,4.5, 
and [Fc/H] — -2.0, -1.0, 0.0 in the temperature range 6750 < 
Tea < 12000 measured in steps of 250 K to T cff = 10 000 
and 500 K thereafter. The parameter space corresponds to 
—0.1 < (B — V)o < 0.3, appropriate for our observations. 
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The resultant spectra were measured in the same way as the 
KSK data. 



